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ABSTRACT 

We present optical photometric and spectroscopic observations of supernova 2013ej. It is one of 
the brightest type II supernovae exploded in a nearby 10 Mpc) galaxy NGC 628. The light curve 
characteristics are similar to type II SNe, but with a relatively shorter 85 day) and steeper (~ 1.7 
mag (100 d) -1 in V) plateau phase. The SN shows a large drop of 2.4 mag in V band brightness 
during plateau to nebular transition. The absolute ultraviolet (UV) light curves are identical to SN 
2012aw, showing a similar UV plateau trend extending up to 85 days. The radioactive 56 Ni mass 
estimated from the tail luminosity is O.O2M 0 which is significantly lower than typical type IIP SNe. 

The characteristics of spectral features and evolution of line velocities indicate that SN 2013ej is a type 
II event. However, light curve characteristics and some spectroscopic features provide strong support 
in classifying it as a type IIL event. A detailed SYNOW modelling of spectra indicates the presence 
of some high velocity components in Ha and H/? profiles, implying possible ejecta-CSM interaction. 

The nebular phase spectrum shows an unusual notch in the Ha emission which may indicate bipolar 
distribution of 56 Ni. Modelling of the bolometric light curve yields a progenitor mass of ~ 14M e and 
a radius of ~ 45OR 0 , with a total explosion energy of ~ 2.3 x 10 51 erg. 

Subject headings: supernovae: general — supernovae: individual: SN 2013ej — galaxies: individual: 
NGC 0628 


1. INTRODUCTION 

Type II supernovae (SNe) originate from massive stars 
with Mzams > 8M 0 (Burrows 2013) which have re¬ 
tained substantial hydrogen in the envelope at the time 
of explosion. They belong to a subclass of core-collapse 
SNe (CCSNe), which collapse under their own gravity at 
the end of the nuclear burning phase, having insufficient 
thermal energy to withstand the collapse. 

The most common subtype among hydrogen rich su¬ 
pernovae is type IIP. At the time of shock breakout al¬ 
most the entire mass of hydrogen is ionized. Type IIP 
SNe have an extended hydrogen envelope, which recom¬ 
bines slowly over a prolonged duration sustaining the 
plateau phase. During this phase the SN light curve 
shows almost constant brightness lasting for 80-100 days. 
At the end of plateau phase the SN experiences a sud- 
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den drop in luminosity, settling onto the slow declining 
radioactive tail, also known as nebular phase, which is 
mainly powered by gamma rays released from the decay 
of 56 Co to 56 Fe, which in turn depends upon the amount 
of 56 Ni synthesized at the time of explosion. 

The plateau slope of SN type II light curve primarily 
depends on the amount of hydrogen present in the ejecta. 
If hydrogen content is high, as in type IIP, the initial en¬ 
ergy deposited from shock and decay of freshly produced 
56 Ni shall be released slowly over a longer period of time. 
On the other hand if hydrogen content is relatively low, 
the light curve will decline fast but with higher peak lu¬ 
minosity. Thus if hydrogen content is low enough, one 
would expect a linear decline in the light curve classifying 
it as type IIL. By the historical classification, type IIL 
(Barbon et al. 1979) shows linear decline in light curve 
over 100 days until it reaches the radioactive tail phase. 
Arcavi et al. (2012) claimed to find type IIP and IIL as 
to distinct group of events which may further indicate 
their distinct class of progenitors. However, recent stud¬ 
ies by Anderson et al. (2014b) and Sanders et al. (2015) 
on large sample of type II SNe do not favor any such bi¬ 
modality in the diversity, rather they found continuum in 
light curve slopes as well as in other physical parameters. 
The continuous distribution of plateau slopes in type II 
events is rather governed by variable amount of hydrogen 
mass left in the envelope at the time of explosion. Based 
on a sample of 11 type IIL events, Faran et al. (2014) 
proposed that any event having decline of 0.5mag in V 
band light curve in first 50 days can be classified as type 
IIL. In light of these recent developments a large num¬ 
ber of type IIP SNe classified earlier may now fall under 
IIL class. Thus many of the past studies collectively on 
samples of type IIP SNe, which we shall be referring in 
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this work may include both IIP as well as IIL. 

Extensive studies have been done to relate observable 
parameters and progenitor properties of IIP SNe (e.g., 
Litvinova & Nadezhin 1985; Hamuy 2003). Stellar evo¬ 
lutionary models suggest that these SNe may originate 
from stars with zero-age-main-sequence mass of 9-25M© 
(e.g., Heger et al. 2003). However, progenitors directly 
recovered for a number of nearby IIP SNe, using the pre- 
SN HST archival images, are found to he within 8 —17M© 
RSG stars (Smartt 2009). Recent X-ray study also infers 
an upper mass limit of < 19M© for type IIP progenitors 
(Dwarkadas 2014), which is in close agreement to that 
obtained from direct detection of progenitors. 

The geometry of the explosion and presence of pre¬ 
existent circumstellar medium (CSM), often associated 
with progenitor mass loss during late stellar evolution¬ 
ary phase, can significantly alter the observables even 
though originating from similar progenitors. There 
are number of recent studies of II SNe, like 2007od 
(Inserra et al. 2011), 2009bw (Inserra et al. 2012) and 
2013by (Valenti et al. 2015) which show signature of such 
CSM interactions during various phases of evolution. 

SN 2013ej is one of the youngest detected type II SN 
which was discovered soon after its explosion. The ear¬ 
liest detection was reported on July 24.125 UTC, 2013 
by C. Feliciano in Bright Supernovae 10 and subsequent 
independent detection on July 24.83 UTC by Lee et al. 
(2013) at V-band magnitude of ^14.0. The last non¬ 
detection was reported on July 23.54 UTC, 2013 by All 
Sky Automated Survey for Supernovae (Shappee et al. 
2013) at a V-band detection limit of > 16.7 mag. There¬ 
fore, we adopt an explosion epoch (Od) of July 23.8 UTC 
(JD = 2456497.3 =b 0.3), which is chosen in between the 
last non-detection and first detection of SN 2013ej. This 
being one of the nearest and brightest events, it provides 
us with an excellent opportunity to study the origin and 
evolution of type II SN. Some of the basic properties of 
SN 2013ej and its host galaxy are listed in Table 1. 

Valenti et al. (2014) presented early observations of 
SN 2013ej and using temperature evolution for the first 
week, they estimated a progenitor radius of 400-600 R©. 
Fraser et al. (2014) used high resolution archival images 
from HST to examine the location of SN 2013ej and iden¬ 
tified the progenitor candidate to be a supergiant of mass 
8 — 15.5M©. Leonard et al. (2013) reported unusually 
high polarization using spectropolarimetric observation 
for the week old SN, as implying substantial asymmetry 
in the scattering atmosphere of ejecta. X-ray emission 
has also been detected by Swift XRT (Margutti et al. 
2013), which may indicate SN 2013ej has experienced 
CSM interaction. 

In this work we present photometric and spectroscopic 
observation of SN 2013ej, and carry out qualitative as 
well as quantitative analysis of the various observables 
through modelling and comparison with other archety¬ 
pal SNe. The paper is organized as follows. In section 2 
we describe photometric and spectroscopic observations 
and data reduction. The estimation of line of sight ex¬ 
tinction is discussed in section 3. In section 4 we ana¬ 
lyze the light curves, compare absolute magnitude light 
curves and color curves. We also derive bolometric lumi¬ 
nosities and estimate nickel mass from the tail luminos- 

10 http://www.rochesterastronomy.org/supernova.html 


TABLE 1 

Relevant parameters for the host galaxy NGC 0628 and 
SN 2013ej. 


Parameters 

Value 

Ref. a 

NGC 0628: 

Alternate name 

M74 

2 

Type 

Sc 

2 

RA (J2000) 

a = 01 h 36 m 41?77 

2 

DEC (J2000) 

S = 15°46'59"8 

2 

Abs. Magnitude 

Mb = —20.72 mag 

2 

Distance 

D = 9.6 ±0.7 Mpc 

1 

Distance modulus 

fi = 29.90 ± 0.16 mag 


Heliocentric Velocity 

c^helio — 658 ± 1 km s —1 

2 

SN 2013ej: 

RA (J2000) 

a = 01 h 36 m 48?16 

3 

DEC (J2000) 

5 = 15°45'31"3 


Galactocentric Location 

1'33"E, 2'15"S 


Date of explosion 

t 0 =23.8 July 2013 (UT) 

1 

Reddening 

(JD 2456497.3 ± 0.3) 

E(B -V) = 0.060 ± 0.001 mag 

1 


(1) This paper; (2) HyperLEDA - http://leda.univ-lyonl.fr; 
(3) Kim et al. (2013) 


ity. Optical spectra are analyzed in section 5, where we 
model and discuss evolution of various spectral features 
and compare velocity profile with other type II SNe. In 
section 7, we model the bolometric light curve of SN 
2013ej and estimate progenitor and explosion parame¬ 
ters. Finally in section 8, we summarize the results of 
this work. 

2. OBSERVATION AND DATA REDUCTION 
2.1. Photometry 

Broadband photometric observations in UB VRI filters 
have been carried out from 2.0m IIA Himalayan Chan¬ 
dra Telescope (HCT) telescope at Hanle and ARIES 1.0m 
Sampurananand (ST) and 1.3m Devasthal Fast Optical 
(DFOT) telescopes at Nainital. Additionally SN 2013ej 
has been also observed with Swift Ultraviolet/optical 
(UVOT) telescope in all six bands. 

Photometric data reductions follows the same pro¬ 
cedure as described in Bose et al. (2013). Images are 
cleaned and processed using standard procedures of 
IRAF software. DAOPHOT routines have been used 
to perform PSF photometry and extracting differential 
light-curves. To standardize the SN field, three Landolt 
standard fields (PG 0231, PG 2231 and SA 92) were ob¬ 
served on October 27, 2013 with 1.0-m ST under good 
photometric night and seeing (typical FWHM ^2".l in 
V band) condition. For atmospheric extinction measure¬ 
ment, PG 2231 and PG 0231 were observed at different 
air masses. The SN field has been also observed in be¬ 
tween standard observations. The standardization coef¬ 
ficients derived are represented in the following transfor¬ 
mation equations, 

u = U + (7.800 =b 0.005) - (0.067 =b 0.009) • (U - B) 
b = B + (5.269 =b 0.007) - (0.060 =b 0.009) • (B - V) 
v = V + (4.677 ± 0.004) - (0.056 ± 0.005) • (B - V) 






Type IIL Supernova 2013ej 


3 


TABLE 2 

Eight local standards in the field of SN 2013ej with corresponding coordinates (a, 6) calibrated magnitudes in UBVRI. 
Errors quoted here, include both photometric and calibration errors. 


Star 

ID 

<AJ2000 

(h m s) 

^ J2000 

(O'") 

V 

(mag) 

B 

(mag) 

V 

(mag) 

R 

(mag) 

/ 

(mag) 

A 

1:36:57.9 

+15:51:19.4 

16.773 ± 0.0325 

16.867 ± 0.0259 

16.297 ± 0.0193 

15.939 ± 0.0163 

15.567 ± 0.0207 

B 

1:36:23.0 

+15:47:45.3 

15.102 ± 0.0289 

15.109 ± 0.0302 

14.580 ± 0.0234 

14.253 ± 0.0183 

13.888 ± 0.0260 

C 

1:36:50.4 

+15:40:01.9 

16.588 ± 0.0318 

16.525 ± 0.0277 

15.798 ± 0.0200 

15.372 ± 0.0157 

14.925 ± 0.0187 

D 

1:36:52.7 

+15:40:39.4 

14.811 ± 0.0318 

14.561 ± 0.0265 

13.817 ± 0.0184 

13.384 ± 0.0167 

12.976 ± 0.0196 

E 

1:37:03.4 

+15:41:39.2 

17.140 ± 0.0386 

17.064 ± 0.0251 

16.407 ± 0.0200 

16.008 ± 0.0160 

15.601 ± 0.0206 

F 

1:37:09.0 

+15:41:20.4 

18.804 ± 0.1537 

17.800 ± 0.0282 

16.769 ± 0.0249 

16.146 ± 0.0167 

15.583 ± 0.0205 

G 

1:36:57.6 

+15:46:22.7 

13.934 ± 0.0272 

13.756 ± 0.0219 

12.991 ± 0.0160 

12.555 ± 0.0161 

12.155 ± 0.0240 

H 

1:37:09.0 

+15:48:00.6 

16.974 ± 0.0434 

16.172 ± 0.0249 

15.175 ± 0.0157 

14.598 ± 0.0170 

14.062 ± 0.0194 



Fig. 1.— SN 2013ej in NGC 0628. The Bl+band composite 
image taken from 104-cm Sampurnanad telescope, covering an area 
of about 13' X13' is shown. Eight local field standards and SN are 
marked in the image. 

r = R+ (4.405 ± 0.005) - (0.038 ± 0.010) • (V - R) 
i = I + (4.821 ± 0.006) - (0.048 =b 0.006) • (V — I) 

where u, 6, t?, r and i are instrumental magnitudes cor¬ 
rected for time, aperture and airmass; 1/, F>, V, R and I 
are standard magnitude. The standard-deviation of the 
difference between the calibrated and the standard mag¬ 
nitudes of the observed Landolt stars are found to be ^ 
0.03 mag in 1/, ~ 0.02 mag in BR and ^ 0.01 mag in VI. 
The transformation coefficients were then used to gener¬ 
ate eight local standard stars in the field of SN 2013ej, 
which are verified to be non-variable and have bright¬ 
ness similar to SN. These stars are identified in Fig.l 
and the calibrated UBVRI magnitudes are listed in Ta¬ 
ble 2. These selected eight local standards were further 
used to standardize the instrumental light curve of the 
SN. One of these stars (star B) is common to that used 
in the study by Richmond (2014), and its BVRI magni¬ 
tudes are found to he within 0.03 mag of our calibrated 
magnitudes. Our calibrated magnitudes for SN 2013ej 
are also found to be consistent within errors to that pre¬ 
sented in earlier studies of the event (Valenti et al. 2014; 
Richmond 2014). The standard photometric magnitudes 
of SN 2013ej are listed in Table 3. 


This supernova was also observed with the Ultra- 
Violet/Optical Telescope (UVOT; Roming et al. 2005) in 
six bands (viz. uvw2, uvm2, uvwl, uvu, uvb, uvv) on 
the Swift spacecraft (Gehrels et al. 2004). The UV pho¬ 
tometry was obtained from the Swift Optical/Ultraviolet 
Supernova Archive 11 (SOUSA; Brown et al. 2014). The 
reduction is based on that of Brown et al. (2009), includ¬ 
ing subtraction of the host galaxy count rates and uses 
the revised UV zeropoints and time-dependent sensitiv¬ 
ity from Breeveld et al. (2011). The UVOT photometry 
is listed in Table. 3. The first month of UVOT photom¬ 
etry was previously presented by Valenti et al. (2014). 

2.2. Spectroscopy 

Spectroscopic observations have been carried out at 10 
phases during 12 to 125d. Out of these, nine epochs of 
low resolution spectra are obtained from Himalaya Faint 
Object Spectrograph and Camera (HFOSC) mounted on 
2.0m HCT. Spectroscopy on the HCT/HFOSC was done 
using a slit width of 1.92 arcsec, and grisms with res¬ 
olution A/AA = 1330 for Gr7 and 2190 for Gr8, and 
bandwidth coverage of 0.38 — 0.64 fim and 0.58 — 0.84 
tim respectively. One high resolution spectrum is ob¬ 
tained from the ARC Echelle Spectrograph (ARCES) 
mounted on 3.5m ARC telescope located at Apache Point 
Observatory (APO). ARCES is a high resolution cross¬ 
dispersion echelle spectrograph, the spectrum is recorded 
in 107 echelle orders covering a wavelength range of A ^ 
0.32-1.00/im, at resolution of R ~ 31500 (Wang et al. 
2003). Summary of spectroscopic observations is given 
in Table. 4. 

Spectroscopic data reduction was done under the IRAF 
environment. Standard reduction procedures are fol¬ 
lowed for bias subtraction and flat fielding. Cosmic ray 
rejections are done using a Laplacian kernel detection 
algorithm for spectra, L.A.Cosmic (van Dokkum 2001). 
One dimensional low resolution spectra were extracted 
using the apall task. Wavelength calibration was done 
using the identify task applied on FeNe and FeAr (for 
HCT) arc spectra taken during observation. Wavelength 
calibration was crosschecked against the [OI] A5577 
sky line in the sky spectrum, and it was found to he 
within 0.3 to 4.5 A of the actual value. Spectra were flux 
calibrated using standard, sensfunc and calibrate 
tasks in IRAF. For flux calibration, spectrophotometric 
standards were used which were observed on the same 
nights as the SN spectra were recorded. All spectra were 
tied to absolute flux scale using the observed flux from 
UBVRI photometry of SN. To perform the tying, indi- 

11 http: //swift, gsfc.nasa.gov/docs/swift / sne/swift _sn.html 
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TABLE 3 

Photometric evolution of SN 2013ej. Errors denote la uncertainty. 

UBVRI photometry 

UT Date JD Phase a U B V R I Tel 6 


(yyyy-mm-dd) 2456000+ (day)(mag)(mag)(mag)(mag)(mag) 


2013-08-04.82 

509.32 

12.02 

12.026 ± 

0.061 

12.633 

± 

0.020 

12.612 

± 

0.013 

12.434 

± 

0.017 

12.349 ± 

0.018 

HCT 

2013-08-31.93 

536.43 

39.13 

14.576 ± 

0.251 

14.208 

± 

0.020 

13.125 

± 

0.011 

12.670 

± 

0.015 

12.436 ± 

0.011 

HCT 

2013-09-29.77 

565.27 

67.97 

16.088 ± 

0.027 

14.991 

± 

0.020 

13.569 

± 

0.012 

13.056 

± 

0.016 

12.750 ± 

0.016 

HCT 

2013-09-30.72 

566.22 

68.92 

16.207 ± 

0.109 

14.956 

± 

0.020 

13.595 

± 

0.008 

12.992 

± 

0.015 

12.741 

± 

0.016 

ST 

2013-10-02.87 

568.37 

71.07 

16.223 ± 

0.028 

15.017 

± 

0.020 

13.640 

± 

0.012 

13.053 

± 

0.015 

12.750 ± 

0.017 

HCT 

2013-10-13.70 

579.20 

81.90 

16.823 ± 

0.054 

15.291 

± 

0.015 

13.864 

± 

0.010 

13.222 

± 

0.017 




ST 

2013-10-15.85 

581.35 

84.05 

17.026 ± 

0.061 

15.365 

± 

0.021 

13.884 

± 

0.012 

13.273 

± 

0.011 

12.978 

± 

0.017 

ST 

2013-10-16.71 

582.21 

84.91 

17.036 ± 

0.089 

15.406 

± 

0.013 

13.939 

± 

0.010 

13.288 

± 

0.018 

12.986 

± 

0.025 

ST 

2013-10-21.73 

587.23 

89.93 

17.292 ± 

0.057 

15.611 

± 

0.017 

14.126 

± 

0.017 

13.446 

± 

0.016 

13.147 

± 

0.018 

ST 

2013-10-24.70 

590.20 

92.90 

17.405 ± 

0.035 

15.743 

± 

0.016 

14.233 

± 

0.021 

13.540 

± 

0.014 

13.253 

± 

0.014 

ST 

2013-10-25.72 

591.22 

93.92 

17.365 ± 

0.023 

15.732 

± 

0.014 

14.340 

± 

0.008 

13.592 

± 

0.011 

13.322 

± 

0.011 

DFOT 

2013-10-26.74 

592.24 

94.94 

17.442 ± 

0.020 

15.795 

± 

0.014 

14.431 

± 

0.007 

13.672 

± 

0.010 

13.384 

± 

0.011 

DFOT 

2013-10-27.76 

593.26 

95.96 

17.515 ± 

0.033 

15.985 

± 

0.022 

14.453 

± 

0.016 

13.750 

± 

0.016 

13.447 

± 

0.021 

ST 

2013-11-09.63 

606.13 

108.83 

18.440 ± 

0.039 

17.611 

± 

0.015 

16.108 

± 

0.012 

15.144 

± 

0.015 

14.783 

± 

0.016 

HCT 

2013-11-11.72 

608.22 

110.92 

18.655 ± 

0.106 

17.725 

± 

0.020 

16.358 

± 

0.012 

15.357 

± 

0.016 

14.978 

± 

0.016 

ST 

2013-11-12.67 

609.17 

111.87 

— 


17.700 

± 

0.021 

16.379 

± 

0.014 

15.358 

± 

0.017 

15.004 

± 

0.014 

ST 

2013-11-14.65 

611.15 

113.85 

— 


17.764 

± 

0.031 

16.405 

± 

0.011 

15.402 

± 

0.010 

15.031 

± 

0.013 

ST 

2013-11-19.69 

616.19 

118.89 

18.515 ± 

0.133 

17.865 

± 

0.023 

16.493 

± 

0.015 

15.480 

± 

0.016 

15.133 

± 

0.019 

ST 

2013-11-23.69 

620.19 

122.89 

19.144 ± 

0.408 

17.945 

± 

0.021 

16.533 

± 

0.009 

15.529 

± 

0.010 

15.203 

± 

0.011 

ST 

2013-11-24.62 

621.12 

123.82 

18.973 ± 

0.128 

17.911 

± 

0.019 

16.552 

± 

0.012 

15.544 

± 

0.015 

15.205 

± 

0.016 

ST 

2013-12-06.72 

633.22 

135.92 

19.292 ± 

0.171 

18.113 

± 

0.028 

16.771 

± 

0.014 

15.719 

± 

0.016 

15.420 

± 

0.017 

ST 

2013-12-08.73 

635.23 

137.93 

19.286 ± 

0.175 

18.139 

± 

0.018 

16.815 

± 

0.017 

15.766 

± 

0.022 

15.486 

± 

0.024 

ST 

2013-12-09.69 

636.19 

138.89 

— 


18.167 

± 

0.022 

16.832 

± 

0.011 

15.779 

± 

0.017 

15.488 

± 

0.017 

ST 

2013-12-10.61 

637.11 

139.81 

— 


18.209 

± 

0.034 

16.863 

± 

0.019 

15.796 

± 

0.020 

15.490 

± 

0.022 

ST 

2013-12-14.74 

641.24 

143.94 

— 


18.015 

± 

0.093 

16.892 

± 

0.034 

15.856 

± 

0.020 

15.597 

± 

0.023 

ST 

2013-12-15.63 

642.13 

144.83 

— 


18.223 

± 

0.041 

16.974 

± 

0.019 

15.914 

± 

0.025 

15.603 

± 

0.026 

ST 

2013-12-16.70 

643.20 

145.90 

— 


18.109 

± 

0.053 

16.943 

± 

0.025 

15.903 

± 

0.019 

15.596 

± 

0.126 

ST 

2013-12-19.61 

646.11 

148.81 

— 


18.249 

± 

0.043 

17.009 

± 

0.015 

15.932 

± 

0.019 

15.661 

± 

0.023 

ST 

2013-12-24.62 

651.12 

153.82 

19.474 ± 

0.061 

18.265 

± 

0.027 

17.138 

± 

0.014 

16.003 

± 

0.015 

15.743 

± 

0.016 

ST 

2013-12-25.66 

652.16 

154.86 

— 


18.321 

± 

0.016 

17.101 

± 

0.010 

16.012 

± 

0.009 

15.722 

± 

0.012 

ST,DFOT 

2013-12-28.62 

655.12 

157.82 

19.368 ± 

0.058 

18.325 

± 

0.019 

17.161 

± 

0.009 

16.041 

± 

0.015 

15.760 

± 

0.016 

DFOT 

2013-12-29.59 

656.09 

158.79 

19.436 ± 

0.060 

18.315 

± 

0.024 

17.180 

± 

0.011 

16.061 

± 

0.010 

15.791 

± 

0.011 

DFOT 

2014-01-19.62 

677.12 

179.82 

— 


18.676 

± 

0.025 

17.458 

± 

0.011 

16.370 

± 

0.014 

16.128 

± 

0.015 

ST 

2014-01-25.62 

683.12 

185.82 

19.703 ± 

0.071 

18.638 

± 

0.013 

17.526 

± 

0.009 

16.424 

± 

0.011 

16.164 

± 

0.012 

DFOT 

2014-01-30.62 

688.12 

190.82 

19.797 ± 

0.596 

18.785 

± 

0.027 

17.602 

± 

0.014 

16.501 

± 

0.013 

16.282 

± 

0.015 

ST 

2014-01-31.58 

689.08 

191.78 

— 


18.787 

± 

0.030 

17.618 

± 

0.019 

16.522 

± 

0.017 

16.273 

± 

0.025 

ST 

2014-02-02.62 

691.12 

193.82 

— 


18.813 

± 

0.035 

17.623 

± 

0.031 

16.546 

± 

0.020 

16.323 

± 

0.024 

ST 

2014-02-17.59 

706.09 

208.79 

— 


19.218 

± 

0.079 

17.814 

± 

0.022 

16.682 

± 

0.012 

16.470 

± 

0.017 

ST 


Swift UVOT photometry 


UT Date 

JD 

Phase a 

uvw 2 

uvm2 

uvw 1 

uvu 

uvb 

uvv 

Tel 6 

(yyyy/mm/dd) 

2456000+ 

(day) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

/Inst 

2013-07-30.98 

504.48 

7.18 

12.369 ± 

0.040 

12.023 ± 0.040 

11.711 ± 0.039 




— 

12.689 ± 0.042 

Swift 

2013-07-31.50 

505.00 

7.70 

12.455 ± 

0.040 

12.097 ± 0.040 

11.755 ± 0.039 




— 

12.614 ± 0.040 

Swift 

2013-07-31.83 

505.33 

8.03 

12.577 ± 

0.035 

12.204 ± 0.033 

11.814 ± 0.032 




— 

— 

Swift 

2013-08-03.06 

507.56 

10.26 

13.044 ± 

0.037 

12.695 ± 0.041 

— 

11.675 ± 

0.029 

12.619 ± 0.029 

— 

Swift 

2013-08-03.18 

507.68 

10.38 

13.056 ± 

0.035 

— 

— 




12.622 ± 0.029 

— 

Swift 

2013-08-04.85 

509.35 

12.05 

13.374 ± 

0.040 

13.155 ± 0.053 

— 

11.812 ± 

0.029 

12.608 ± 0.029 

— 

Swift 

2013-08-04.98 

509.48 

12.18 

13.385 ± 

0.037 

— 

— 




— 

— 

Swift 

2013-08-07.24 

511.74 

14.44 

13.907 ± 

0.041 

— 

12.948 ± 0.042 




— 

— 

Swift 

2013-08-07.55 

512.05 

14.75 

13.968 ± 

0.050 

— 

— 




— 

— 

Swift 

2013-08-08.02 

512.52 

15.22 

14.039 ± 

0.052 

14.058 ± 0.070 

13.131 ± 0.038 

12.185 ± 

0.031 

12.749 ± 0.029 

12.477 ± 0.030 

Swift 

2013-08-08.22 

512.72 

15.42 

14.126 ± 

0.045 

— 

— 

12.266 ± 

0.029 

— 

— 

Swift 

2013-08-09.25 

513.75 

16.45 

14.387 ± 

0.055 

14.305 ± 0.112 

13.379 ± 0.041 

12.333 ± 

0.029 

12.906 ± 0.029 

12.535 ± 0.031 

Swift 

2013-08-09.31 

513.81 

16.51 

— 


14.406 ± 0.065 

— 




— 

— 

Swift 

2013-08-11.78 

516.28 

18.98 

15.210 ± 

0.118 

15.114 ± 0.109 

13.907 ± 0.052 

12.659 

± 

0.029 

12.983 ± 0.029 

12.581 dh 0.031 

Swift 

2013-08-13.85 

518.35 

21.05 

15.652 ± 

0.082 

15.964 ± 0.068 

14.446 ± 0.059 

12.982 

± 

0.031 

13.109 ± 0.029 

12.599 dh 0.032 

Swift 

2013-08-15.00 

520.50 

23.20 

16.209 ± 

0.090 

— 

14.905 ± 0.069 

13.308 

± 

0.033 

13.221 ± 0.030 

12.573 ± 0.032 

Swift 

2013-08-17.65 

522.15 

24.85 

16.588 ± 

0.098 

17.109 ± 0.195 

15.201 ± 0.072 

13.602 

± 

0.035 

13.293 ± 0.030 

12.656 ± 0.032 

Swift 

2013-08-19.73 

524.23 

26.93 

16.824 ± 

0.105 

17.554 ± 0.221 

15.493 ± 0.076 

13.964 

± 

0.039 

13.476 ± 0.030 

12.692 ± 0.032 

Swift 

2013-08-22.54 

527.04 

29.74 

17.245 ± 

0.120 

18.047 ± 0.250 

15.890 ± 0.075 

14.338 

± 

0.045 

13.663 ± 0.032 

12.816 dh 0.033 

Swift 

2013-08-23.14 

527.64 

30.34 

17.170 ± 

0.117 

— 

15.866 ± 0.083 

14.366 

± 

0.045 

13.627 ± 0.031 

12.892 ± 0.034 

Swift 

2013-08-27.74 

532.24 

34.94 

17.746 ± 

0.146 

18.569 ± 0.214 

16.356 ± 0.095 

14.844 

± 

0.058 

13.915 ± 0.033 

12.965 dh 0.034 

Swift 

2013-09-06.16 

541.66 

44.36 

18.133 ± 

0.124 

19.137 ± 0.190 

16.793 ± 0.084 

15.573 

± 

0.067 

— 

— 

Swift 

2013-09-06.41 

541.91 

44.61 

— 


— 

— 

15.674 

± 

0.087 

14.367 ± 0.036 

13.231 ± 0.035 

Swift 

2013-09-16.71 

552.21 

54.91 

18.687 ± 

0.158 

19.486 ± 0.236 

17.292 ± 0.096 

16.229 

± 

0.090 

14.750 ± 0.039 

13.470 ± 0.038 

Swift 

2013-09-26.45 

561.95 

64.65 

18.793 ± 

0.166 

— 

17.562 ± 0.123 

16.585 

± 

0.128 

14.922 ± 0.042 

13.604 ± 0.039 

Swift 

2013-10-06.88 

572.38 

75.08 

19.241 ± 

0.231 

19.883 ± 0.333 

17.919 ± 0.133 

17.055 

± 

0.094 

— 

— 

Swift 

2013-10-16.77 

582.27 

84.97 

19.294 ± 

0.247 

— 

18.127 ± 0.170 

17.286 

± 

0.164 

15.464 ± 0.055 

14.029 ± 0.045 

Swift 

2013-10-26.95 

592.45 

95.15 

— 


— 

18.248 ± 0.190 

17.514 

± 

0.171 

— 

— 

Swift 

2013-11-06.16 

602.66 

105.36 

— 


— 

— 

18.774 

± 

0.304 

— 

— 

Swift 

2013-11-13.21 

609.71 

112.41 

— 


— 

19.523 ± 0.351 

19.058 

± 

0.306 

— 

— 

Swift 

2013-11-13.68 

610.18 

112.88 

— 


— 

— 

18.816 

± 

0.263 

17.974 ± 0.210 

16.512 ± 0.148 

Swift 

2013-11-20.43 

616.93 

119.63 

— 


— 

— 

18.977 

± 

0.214 

17.889 ± 0.090 

16.718 dh 0.080 

Swift 

2013-11-25.40 

621.90 

124.60 

— 


— 

— 

19.162 

± 

0.323 

— 

— 

Swift 

2013-11-30.43 

626.93 

129.63 

— 


— 

19.726 ± 0.313 

19.342 

± 

0.280 

18.155 ± 0.101 

16.834 ± 0.082 

Swift 

2013-12-09.75 

636.25 

138.95 

— 


— 

19.807 ± 0.327 

19.343 

± 

0.274 

18.196 ± 0.102 

16.928 ± 0.085 

Swift 


a with reference to the explosion epoch JD 2456497.30 

6 ST : 104-cm Sampurnanand Telescope, ARIES, India; DFOT : 130-cm Devasthal fast optical telescope, ARIES, India; HCT: 2m Himalyan Chandra 
Telescope, Hanle, India; Swift: Swift UVOT 

Note: Data observed within 5 Hrs, are represented under single epoch observation. 
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TABLE 4 

Summary of spectroscopic observations of SN 2013ej. The spectral observations are made at 10 phases during 12d to 125d. 


UT Date 
(yy/mm/dd.dd) 

JD 

2456000+ 

Phase® 

(days) 

Telescope 0 

Range b 

fi m 

Exposure 

(s) 

2013-08-04.86 

509.36 

12.1 

HUT 

0.38-0.68 

900 

2013-08-27.76 

532.26 

35.0 

HCT 

0.38-0.68 

1200 




HCT 

0.58-0.84 

1200 

2013-09-03.90 

539.40 

42.1 

HCT 

0.38-0.68 

1500 




HCT 

0.58-0.84 

1500 

2013-09-29.78 

565.28 

68.0 

HCT 

0.38-0.68 

1800 




HCT 

0.58-0.84 

2400 

2013-10-02.89 

568.39 

71.1 

HCT 

0.38-0.68 

1500 




HCT 

0.58-0.84 

1500 

2013-10-11.28 

576.78 

79.5 

APO 

0.32-1.00 

1200 

2013-10-27.87 

593.37 

96.1 

HCT 

0.38-0.68 

2400 

2013-10-28.79 

594.29 

97.0 

HCT 

0.58-0.84 

2400 

2013-11-09.65 

606.15 

108.9 

HCT 

0.38-0.68 

2100 




HCT 

0.58-0.84 

3900 

2013-11-25.75 

622.25 

125.0 

HCT 

0.38-0.68 

2400 




HCT 

0.58-0.84 

2400 


a With reference to the adopted explosion time JD 2456497.30 
b For transmission >50% 

c HCT : HFOSC on 2 m Himalyan Chandra Telescope, India; APO : Echelle spectrograph on 3.5 m ARC telescope at Apache Point 
Observatory, U.S. 
d At 0.6 ^m 
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1 _1 

79.5d 

Na ID : M.W Na ID: NGC 0628 


5880 5885 5890 5895 5900 5905 5910 5915 

Observed wavelength (A) 


Fig. 2.— Echelle spectra at 79.5d showing the Nai D doublet 
for Milky-way while no impression for NGC 0628 is detected. 


vidual spectrum is multiplied by a wavelength dependent 
polynomial, which is convolved with UB VRI filters and 
then the polynomial is tuned to match the convolved flux 
with observations. The one dimensional calibrated spec¬ 
tra were corrected for heliocentric velocity of host galaxy 
(658 kms -1 ; Table 1) using DOPCOR task. 

3. DISTANCE AND EXTINCTION 

We adopt a distance of 9.57 =b 0.70 Mpc which is a 
mean value of four different distance estimation tech¬ 
niques used for NGC 0628, viz., 9.91 Mpc apply¬ 
ing Standard Candle Method (scm) to SN 2003gd by 
Olivares et al. (2010); 10.19 Mpc using the Tully-Fisher 
method (HyperLeda 12 ); 9.59 Mpc using brightest su¬ 
pergiant distance estimate by Hendry et al. (2005); and 
planetary nebula luminosity function distance 8.59 Mpc 
(Herrmann et al. 2008). Although for each of these 
methods number of distance estimates exists in liter¬ 
ature, we tried to select only most recent estimates. 
Richmond (2014) estimated a distance of 9.1 ± 0.4 Mpc 
by applying Expanding Photosphere Method (EPM) to 
SN 2013ej, which we find consistent to that we adopted. 

One of the most reliable and well accepted method 
for SNe line-of-sight reddening estimation is using the 
Nai D absorption feature. The equivalent width 
(EW) of Nai D doublet (AA 5890, 5896) is found to 
be correlated with the reddening, estimated from the 
tail color curves of type la SNe (Barbon et al. 1990; 
Turatto et al. 2003). However, Poznanski et al. (2011) 
suggested that although Nai D EW is weakly correlated 
with E(B — V), the EWs estimated from low resolution 
spectra is a bad estimator of E(B — V). Poznanski et al. 
(2012) used a larger sample of data and presented a more 
precise and rather different functional form of the corre¬ 
lation than that was derived earlier. Our high resolution 
echelle spectra at 79.5d provided an excellent opportu¬ 
nity to investigate the line-of-sight extinction. 

The resolved Nai D doublet for Milky-way is clearly 
visible in the high-resolution spectra (recorded on 79.5d) 
as shown in Fig. 2. Whereas no impression of Nai D 
for NGC 0628 is detected at the expected redshifted po¬ 
sition relative to Milky-way. This indicates that the 
reddening due to host is negligible, only Galactic red¬ 
dening will contribute to the total line of sight extinc¬ 
tion. A similar conclusion has also been inferred by 

12 http://leda.univ-lyonl.fr/ 



Fig. 3. — The photometric light curves in Johnson-Cousins UB- 
VRI and Swift UVOT bands. The light curves are vertically shifted 
for clarity. The line joining the data points of light curves is for 
visualization purpose only. 


Valenti et al. (2014) from their high resolution spectra 
obtained at 3Id. Thus, we adopt a total E(B — V) = 
0.060 ±0.001 mag, which is entirely due to Galactic red¬ 
dening (Schlafly & Finkbeiner 2011) and assuming total- 
to-selective extinction at V band as Ry = 3.1, it trans¬ 
lates into Ay = 0.185 ± 0.004 mag. 

4. LIGHT CURVE 

4.1. Light curve evolution and comparison 

The optical light curves of SN 2013ej in UBVRI and six 
UVOT bands are shown in Fig. 3. UBVRI photometric 
observations were done at 38 phases during 12 to 209d 
(from plateau to nebular phase). The duration of plateau 
phase is sparsely covered, while denser follow-up initiated 
after 68d. The plateau phase lasted until ~ 85d with an 
average decline rate of 6.60, 3.57, 1.74, 1.07 and 0.74 
mag (100 d) -1 in UBVRI bands respectively. Since 95d, 
the light curve declines very fast until 115d, after which 
it settles to a relatively slow declining nebular phase. 
During this phase the decline rates for UBVRI bands are 
0.98, 1.22, 1.53, 1.42 and 1.55 mag (100 d) -1 respectively. 

SN 2013ej has been also observed by Swift UVOT at 
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35 phases during 7 to 139d. The UVOT UV band light 
curves declines steeply during the first 30d at a rate of 
0.182, 0.213, 0.262 mag d -1 in uvwl, uvw2 and uvm2 
bands respectively, thereafter settling into a slow declin¬ 
ing phase until it reaches the end of plateau. 

SN 2013ej experience a steeper plateau decline than 
that observed for SN 1999em (Leonard et al. 2002c), SN 
1999gi (Leonard et al. 2002b), SN 2012aw (Bose et al. 
2013) and SN 2013ab (Bose et al. 2015). For example, 
SN 2012aw plateau declines at a rate of 5.60, 1.74, 0.55 
mag (100 d) -1 in UBV- bands, similarly for SN 2013ab 
decline rates in UBVRI are 7.60, 2.72, 0.92, 0.59 and 
0.30 mag (100 d) -1 and 0.169, 0.236, 0.257 mag d -1 in 
UVOT uvwl, uvw2 and uvm2 bands (during first 30d). 

The absolute V-band (My) light curve of SN 2013ej is 
plotted in Fig. 4 and is compared with other well studied 
type II SNe (after correcting for extinction and distance). 
In Table 5 we list the plateau slope of all compared type 
II events. The comparison shows that the decline rate 
of SN 2013ej during this phase is highest (1.74 mag (100 
d) -1 ) among most other SNe, except three type IIL SNe 
1980K, 2000dc and 2013by, where SN 1980 is among the 
very first observed prototypical type IIL event. The early 
plateau (< 40d) light curve of SN 2013ej is identical 
to SN 2009bw. However, unlike most other IIP SNe, 
e.g. 2009bw and 2013ab, which becomes flatter during 
late plateau, SN 2013ej continues to decline almost at a 
steady rate until the end of plateau (~ 85d). The mid¬ 
plateau My = —14.7 mag for SN 2013ej, which places it 
in the class of normal luminous type II events. SN 2013ej 
is comparable with fast declining and short plateau SNe 
in the sample of Anderson et al. (2014b). Following the 
plateau phase, V-band light drops very fast to reach slow 
declining nebular phase (1.53 mag (100 d) _1 ), which is 
powered by the radioactive decay of 56 Co to 56 Fe. The 
fall of My during the plateau nebular transition is ~ 2.4 
mag, which is on the higher side of the compared events. 
The closest comparison is SNe 2009bw and 2012A which 
exhibits a drop of ^2.4 mag and ^2.5 mag respectively. 
This also indicates low amount of 56 Ni mass synthesized 
during the explosion which we shall further discuss in the 
next section. 

Swift UVOT absolute magnitude light curves of SN 
2013ej are shown in Fig. 5 and compared with other well 
observed type II SNe. The sample is selected in such a 
way that SNe have at least a month of observations. Most 
SNe are not followed for more than a month by Swift , 
mainly because of the large distances or high extinction 
values. However, both these factors work in favor of SN 
2013ej making it possible to have about four months of 
observations. Moreover, the location of the SN being in 
the outskirt of a spiral arm of NGC 0628, the background 
flux contamination is also negligible. The comparison 
shows that the SN 2013ej UV light curves are identical to 
SN 2012aw. SN 2013ej also shows a similar UV plateau 
trend as observed in SN 2012aw (Bayless et al. 2013), 
which is although expected but rarely detected for IIP/L 
SNe. 

Broadband color provides important information to 
study the temporal evolution of SN envelope. In Fig. 6, 
we plot the intrinsic colors U-B, B-V, V-R and V-I for 
SN 2013ej and compare its evolution with type Il-pec 
SN 1987A, and type IIP SNe 1999em, 2004et, 2012aw 


TABLE 5 

Parameters estimated from V band light cruve 
SN Name Plateau slope® Transition drop 6 Transition time c 


mag (100 d) 1 mag days 


SN1980K 

3.63 

± 

0.04 

2.0T0.04 

37 

± 

5 

SN2000dc 

2.56 

± 

0.06* 

- 


- 


SN2013by 

2.01 

± 

0.02 

2.2T0.03 

19 

± 

5 

SN 2013ej 

1.74 

± 

0.08 

2.4T0.02 

21 

± 

3 

SN2003hn 

1.41 

± 

0.04 

2.0T0.04 

19 

4 = 

4 

SN2012A 

1.12 

± 

0.03 

2.5T0.02 

23 

=b 

4 

SN2009bw 

0.93 

4= 

0.04 

2.4T0.03 

14 

± 

3 

SN2004et 

0.73 

± 

0.02 

2.14=0.04 

27 

± 

6 

SN2013ab 

0.54 

± 

0.02 

1.74=0.02 

25 

4 = 

2 

SN2012aw 

0.51 

4= 

0.02 

- 


- 


SN1999gi 

0.47 

± 

0.02 

2.0T0.02 

29 

4 = 

3 

SN2005cs 

0.44 

± 

0.03 

4.0T0.03 

24 

4 = 

3 

SN2009N 

0.36 

± 

0.03 

2.0T0.04 

26 

± 

3 

SN1999em 

0.31 

± 

0.02 

1.94=0.02 

28 

=b 

4 


Note: Objects are sorted in order of plateau slope. 
a Plateau slope during the linear decline phase, starting after first 
minima until plateau end. 

b Drop in magnitude during the plateau to nebular transition. 
c Duration of plateau to nebular transition. 

* Slope is calculated up to the available range of data, as plateau 
end is not observed. 

and 2013ab. All the colors show generic signature of fast 
cooling ejecta until the end of plateau (~ llOd). With 
the start of the nebular phase it continues to cool at a 
much slower rate in V-I and V-R colors, whereas U-V 
and B-V shows a bluer trend. This is because, as the 
SN enters the nebular phase, the ejecta become depleted 
of free electrons, thereby making the envelope optically 
thin, and so unable to thermalize the photons from ra¬ 
dioactive decay of 56 Co to 56 Fe. 

4.2. Bolometric light-curve 

We compute the pseudo-bolometric luminosities fol¬ 
lowing the method described in Bose et al. (2013); which 
include SED integration over the semi-deconvolved pho¬ 
tometric fluxes after correcting for extinction and dis¬ 
tance. Supernova bolometric luminosities during early 
phases (< 30d) are dominated by ultraviolet fluxes, 
while after mid-plateau (^ 50d) UV contribution be¬ 
comes insignificant as compared to optical counterpart 
(e.g., as seen in SNe 2012aw, 2013ab; Bose et al. 2013, 
2015). Similarly, during late phases > lOOd NIR be¬ 
comes dominant over optical fluxes. However, during 
most of the light curve evolution, optical fluxes still 
provide significant contribution. We compute pseudo- 
bolometric luminosities in the wavelength range of U 
to / band (3335-8750A). We also computed UV-optical 
pseudo-bolometric light curve with wavelength starting 
from uvw2 band (wavelength range of 1606-8750A). The 
UV contribution enhances the luminosity significantly 
during early phases, whereas it is almost negligible af¬ 
ter mid-plateau. 

In Fig. 7, we plot pseudo bolometric light curve for 
SN 2013ej and compare it with other SNe light curves 
computed using the same technique. We also include 
UV-optical bolometric light curve for SNe 2012aw and 
2013ab along with SN 2013ej for comparison. Al¬ 
though the UV-optical light curve is initially brighter 
than the optical light curve, they completely coin¬ 
cide by the end of plateau phase (85d). It is evi¬ 
dent from the comparison that SN 2013ej experienced 
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Fig. 4.— The My light curve of SN 2013ej is compared with other type II SNe. The exponential decline of the tail light curve following 
the radioactive decay law for 56 Co—> 56 Fe is shown with a dashed line. On the bottom left side, pair of dotted lines in each gray and 
green colors represent the slope range for type IIP and IIL SNe templates as given by Faran et al. (2014). The adopted explosion time in 
JD-2400000, distance in Mpc, E(B — V) in mag and the reference for observed V-band magnitude, respectively, are : SN 1980K - 44540.5, 

5.5, 0.30; Barbon et al. (1982), NED database; SN 1987A - 46849.8, 0.05, 0.16; Hamuy & Suntzeff (1990); SN 1999em - 51475.6, 11.7, 
0.10; Leonard et al. (2002a); Elmhamdi et al. (2003); SN 1999gi - 51522.3, 13.0, 0.21; Leonard et al. (2002b); SN 2000dc - 51762.4, 49.0, 
0.07; Faran et al. (2014), NED database; SN 2003hn — 52866.5, 17.0, 0.19; Krisciunas et al. (2009); Anderson et al. (2014b); SN 2004et - 

53270.5, 5.4, 0.41; Sahu et al. (2006); SN 2005cs - 53549.0, 7.8, 0.11; Pastorello et al. (2009); SN 2009N - 54848.1, 21.6, 0.13; Takats et al. 
(2014); SN 2009bw - 54916.5, 20.2, 0.31; Inserra et al. (2012); SN 2012A - 55933.5, 9.8, 0.04; Tomasella et al. (2013); SN 2012aw - 56002.6, 
9.9, 0.07; Bose et al. (2013); SN 2013ab - 56340.0, 24.0, 0.04; Bose et al. (2015); SN 2013by - 56404.0, 14.8, 0.19; Valenti et al. (2015). 
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a steep decline during the plateau phase, but with a 
much shorter duration. This is consistent with the 
anti-correlation observed between plateau slope and du¬ 
ration for type II SNe (Blinnikov & Bartunov 1993; 
Anderson et al. 2014b). The UV-optical bolometric light 
decreases by 0.83 dex during plateau phase (from 12 to 
85d), followed by an even faster drop by 0.76 dex in a 
short duration of 21 days (from 90 to 11 Id). There¬ 
after, the SN settles in a slow declining nebular phase. 
The tail luminosities are significantly lower than other 
normal luminosity IIP events, e.g., SN 2013ej luminosi¬ 
ties are lower by ^ 0.5 dex (at 200d) than that of type 
II SNe 1987A, 1999em, 2004et and 2012aw, but higher 
than subhmhnous events like SN 2005cs. Another notice¬ 
able dissimilarity of the tail light curve is its high decline 
rate. SN 2013ej tail luminosity declines at a rate of 0.55 
dex 100 d -1 , which is much higher than that expected 
from radioactive decay of 56 Co to 56 Fe. This is possibly 
because of inefficient gamma-ray trapping in the ejecta, 
and thus incomplete thermalization of the photons. We 
shall further explore this in §7 in context of modeling the 


light curve. 

4.3. Mass of nickel 

During the explosive nucleosynthesis of silicon and oxy¬ 
gen, at the time of shock-breakout in CCSNe, radioac¬ 
tive 56 Ni is produced. The nebular-phase light-curve is 
mainly powered by the radioactive decay of 56 Ni to 56 Co 
and 56 Co to 56 Fe with half-life times of 6.Id and 77.Id 
respectively emitting y-rays and positrons. Thus the tail 
luminosity will be proportional to the amount of radioac¬ 
tive 56 Ni synthesized at the time of explosion. We deter¬ 
mine the mass of 56 Ni using following two methods. 

For SN 1987A, one of the most well studied and well ob¬ 
served event, the mass of 56 Ni produced in the explosion 
has been estimated quite accurately, to be 0.075 db 0.005 
M 0 (Arnett 1996). By comparing the tail luminosities of 
SN 2013ej and SN 1987A at similar phases, it is possible 
to estimate the 56 Ni mass for SN 2013ej. In principle 
true bolometric luminosities (including UV, optical and 
IR) are to be used for this purpose, which are available 
for SN 1987A, whereas for SN 2013ej we have only UV 





















Type IIL Supernova 2013ej 


9 



0 50 100 0 50 100 

Phase (Days) 


Fig. 5.— Comparison of the Swift UVOT UV absolute light 
curves of SN 2013ej, with other well observed II SNe from 
UVOT. For the compared SNe, references for UVOT data, ex¬ 
tinction and distance are: SN 2005cs - Brown et al. (2009); 
Pastorello et al. (2009), SN 2006at - Brown et al. (2009); Dis¬ 
tance 65 Mpc; E(B — V) = 0.031 mag (only Galactic reddening 
Schlafly & Finkbeiner 2011), SN 2006bp - Dessart et al. (2008), 
SN 2012aw - Bayless et al. (2013); Bose et al. (2013), SN 2013ab - 
Bose et al. (2015), SN 2013by - Brown et al. (2014); Valenti et al. 
(2015). Some late data points for SN 2013ab with large errors has 
been omitted from the plot. 



Phase (Days) 

Fig. 6. — The intrinsic colors evolution of SN 2013ej is com¬ 
pared with other well-studied type IIP SNe 1987A, 1999em, 2004et, 
2012aw and 2013ab. The reference for the data is same as in Fig. 4. 
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Fig. 7. — The UBVRI bolometric light-curve of SN 2013ej is 
compared with other well studied super novae. Light curves with 
added UVOT UV contributions are also shown for SNe 2013ej, 
2013ab and 2012aw (labeled as UVO). The adopted values of dis¬ 
tances, reddening and explosion time are same as in Fig. 4. The 
exponential decline of the tail light curve following the radioactive 
decay law is shown with a dashed line. 


and optical observations. Thus, in order to have unifor¬ 
mity in comparison, we used only the UB VRI bolometric 
luminosities for both SNe and computed using the same 
method and wavelength range. We estimate the tail UB¬ 
VRI luminosity at 175d, by making a linear fit over 155 
to 195d, to be 2.90 =b 0.43 x 10 40 erg s -1 . Likewise, SN 
1987A luminosity is estimated to be 9.60 =b 0.06 x 10 40 
erg s -1 at similar phase. Thus, the ratio of SN 2013ej to 
SN 1987A luminosity is 0.302 ±0.044, which corresponds 
to a 56 Ni mass of 0.023 ± O.OO3M 0 for SN 2013ej. 

Assuming the y-photons emitted from radioactive de¬ 
cay of 56 Co thermalize the ejecta, 56 Ni mass can be inde¬ 
pendently estimated from the tail luminosity as described 
by Hamuy (2003). 


Mni = 7.866 x 10 44 x L t exp 


(U ~ * q )/(1 ± z) - 6.1 

111.26 


M 


where to is the explosion time, 6.Id is the half-life time 
of 56 Ni and 111.26d is the e-folding time of the 56 Co 
decay. We compute tail luminosity L t at 6 epochs 
within 153 to 185d from the V band data corrected 
for distance, extinction and bolometric correction fac¬ 
tor of 0.26 at 0.06 mag during nebular phase (Hamuy 
2003). The weighted mean value of L t is found to be 
5.45 ± 0.35 x 10 40 ergs 1 corresponding to mean phase 
of 170d. This tail luminosity corresponds to a value of 
M m = 0.019 ±0.002M o . 

We take the weighted mean of the estimated values 
from above two methods, and adopt a 56 Ni mass of 
0.020 ± 0.002M o for SN 2013ej. 

Hamuy (2003) found a strong correlation between the 
56 Ni mass and the mid plateau (at 50d) V band abso¬ 
lute magnitude for type II SNe and this correlation was 
further confirmed by Spiro et al. (2014) specifically for 
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Fig. 8. — The plot of absolute V band magnitude at 50 day versus 
56 Ni mass for 34 type II SNe. Data taken from Hamuy (2003) and 
Spiro et al. (2014). The position of SN 2013ej on the correlation is 
shown with filled red circle. 

low luminous events. Fig. 8 shows the correlation of mid 
plateau My versus 56 Ni mass for 34 events, including SN 
2013ej. The SN lies within the scatter relation, but to¬ 
wards the lower mass range of 56 Ni than where most of 
the events cluster around (top right). 

5. OPTICAL SPECTRA 

5.1. Key spectral features 

The spectroscopic evolution of SN 2013ej is presented 
in Fig. 9. Preliminary identifications of spectral features 
has been done as per previously studied type IIP SNe 
(e.g., Leonard et al. 2002a; Bose et al. 2013). The spec¬ 
trum at 12d shows broad Ha, H/3 and Hel features on 
top of a hot blue continuum. The 35d spectrum shows a 
relatively flat continuum with well developed features of 
Ha, H/3, Fe II along with blends of other heavier species 
Till and Ball . Hel line is no longer detectable, in¬ 
stead Nal D features start to appear at similar location. 
The spectra from 35 to 80d represent the cooler photo- 
spheric phase, where the photosphere starts to penetrate 
deeper layers rich in heavier elements like Fe II and Sc II 
. During these phases we see the emergence and develop¬ 
ment of various other heavy atomic lines and their blends 
like Till , Ball , Nal D and Call . Fig. 10 shows the 
comparison of three plateau phase spectra, viz. 12, 35 
and 68d with other well studied type IIP SNe at similar 
epochs. The comparison shows the spectra of SN 2013ej 
is broadly identical to others in terms of observable line 
features and their evolution. A notable feature during 
early spectrum (12d) is the dip on the bluer wing of Ha 
profiles near 6170 A which can be attributed as the Sill 
feature. Leonard et al. (2013) also identified this feature 
at ~ 9d spectra of SN 2013ej however, due to unlikeli¬ 
ness of such a strong Sill feature at such early epochs, 
a possiblity of non-standard red supergiant envelope or 
CSM interaction was suggested. However, such dips are 



detectable in 35 and 42d spectra, which we identify as 
Sill feature in SYNOW modeling. 

The spectra at 96 and 97d represents the plateau- 
nebular transition phase. Thereafter, spectra at 109 and 
125d represents the nebular phase, where the ejecta has 
become optically thin. These spectra shows the emer¬ 
gence of some emission features from forbidden lines of [ 
Oi ] AA 6300, 6364 and [ Can ] AA 7291, 73241, as well 
as previously evolved permitted lines of Hi , and the 
Nal A5893 doublet (see Fig. 11). 

Gutierrez et al. (2014) found correlations between Ha 
absorption to emission strengths and light curve param¬ 
eters, i.e. plateau slope and duration of optically thick 
phase. Following their selection criteria for choosing 
phase of SN spectra, i.e. ten days after start of recombi¬ 
nation, we selected 42d spectrum as the closet available 
phase to the criteria. The Ha absorption to emission 
ratio of equivalent widths for SN 2013ej is found to be 
0.23 zb 0.02, the optically thick phase is ^ 85d and B- 
band late plateau (40 to 85d) slope is ~ 0.27 mag (100 
d) -1 . The correlation for optically thick phase dura¬ 
tion is found to follow that presented by Gutierrez et al. 
(2014). For the plateau slope, the correlation also hold 
true, but here SN 2013ej lies in the border line position 
of the scattered relation. However, it may be noted that 
Ha profiles are possibly contaminated by high velocity 
features as we describe in next sections, which may re¬ 
sult in deviation from correlation. 

5.2. SYNOW modelling of spectra 

SN 2013ej spectra has been modeled with SYNOW 13 
(Fisher et al. 1997, 1999; Branch et al. 2002) for line 
identification and its velocity estimation. SYNOW is 
a highly parametrized spectrum synthesis code which 
employs the Sobolev approximation to simplify radia¬ 
tion transfer equations assuming a spherically symmet¬ 
ric supernova expanding homologously. The strength 
of the SYNOW code is its capability to reproduce P- 
Cygni profiles simultaneously in synthetic spectra for a 
given set of atomic species and ionization states. The 
applicability of SYNOW is well tested in various core¬ 
collapse SNe studies (e.g. Inserra et al. 2012; Bose et al. 
2013; Milisavljevic et al. 2013; Bose & Kumar 2014; 
Takats et al. 2014; Marion et al. 2014) for velocity es¬ 
timation and analysis of spectral lines. 

To model the spectra we tried various optical depth 
profiles (viz. gaussian, exponential and power law) with 
no significant difference among them, however we find 
exponential profile (r oc exp[—v/v e ]) marginally better 
suited to match the absorption trough of observed spec¬ 
tra, where v e the e-folding velocity, is a fitted parameter. 
While modeling spectra, Hi lines are always dealt as 
detached scenario. This implies the velocity of hydrogen 
layer is significantly higher and is thus detached from 
photospheric layer, close to which most heavier atomic 
lines form, as assumed in SYNOW code. As a conse¬ 
quence to this, the Ha lines in synthetic spectrum, which 
are highly detached, has flat topped emissions with blue 
shifted absorption counter parts. 

SN 2013ej spectra are dereddened and approximate 
blackbody temperature is supplied in the model to match 

13 https://c3.1bl.gov/es/#id22 
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Fig. 9.— The redshift corrected spectra of SN 2013ej are plotted for 10 phases during 12d to 125d. The prominent P-Cygni profiles of 
hydrogen (Ha, H/3, Hy) and helium ( Hel A5876) are marked. The telluric absorption features of O 2 are marked with ®, symbol. Portion 
of spectra in extreme blue or red ends have low SNR. Individual spectra with with overall low SNR has been binned for better visualization. 


the spectral continuum. For early spectrum (12d), lo¬ 
cal thermodynamic equilibrium (LTE) assumption holds 
good and thus SYNOW could fit the continuum well, 
whereas at later epochs it fails to fit properly. The set 
of atomic species incorporated to generate the synthetic 
model spectrum are Hi , Hel ; Fell ; Till ; Sell ; 
Call; Ball; Nal and Sill. The photospheric velocity 
u p h is optimized to simultaneously fit the Fe II (AA 4924, 
5018, 5169) P-Cygni profiles and Hi lines are treated as 
detached. The optical depths and optical depth profile 
parameters, e-folding velocity are varied for individual 
species to fit respective line profiles. In Fig. 12 we show 
the model fit of 71d spectrum. Most of the observable 
spectral features are reproduced well and are identified 
in the figure. 

Similarly all spectra during 12 to 97d are modeled with 
SYNOW. The model fits for Fell (AA 4924, 5018, 5169), 
H/3 and Ha spectral sections are shown in Fig. 13. The 
atomic species which are important to model these fea¬ 
tures are Hi , Fell , Ball , Till , Sell and Nal . 
In addition to these Si II is also used to model the dips 
in the blue wing of Ha P-Cygni during 12 to 42d. While 
modeling the Ha and H/3 profiles, SYNOW was unable 
to properly fit the broad and extended P-Cygni absorp¬ 
tion troughs with single regular component. In order to 
fit these extended troughs, we invoke high-velocity (HV) 
component of Hi . Although no separate dip is seen, 


possibly due to low spectral resolution and overlapping 
of broad P-Cygni profiles, the HV component can well 
reproduce the observed features in synthetic model spec¬ 
trum. The implication and interpretation of these HV 
components are further discussed in §5.4. The SYNOW- 
derived velocities for Fell , Hcq H/3 lines and corre¬ 
sponding HV components are listed in Table 6. The 
nebular spectra during 109 to 125d have not been mod¬ 
eled primarily due to limitations of the LTE assumption 
of SYNOW, and also because nebular phase spectra are 
dominated by emission lines rather than P-Cygni pro¬ 
files. 

5.3. Evolution of spectral lines 

Investigation of the spectral evolution sheds light on 
various important aspects of the SN, like interaction of 
ejecta with the circumstellar material, geometrical distri¬ 
bution of expanding shell of ejecta and formation of dust 
during late time. SN spectra are dominated by P-Cygni 
profiles which are direct indicators of expansion veloci¬ 
ties and they evolve with the velocity of photosphere. As 
ejecta expands and opacity decreases allowing photons 
to escape from deeper layers rich in heavier elements, we 
are able to see emergence and growth of various spectral 
lines. 

To illustrate the evolution of Ha line, in Fig. 14 
partial region of spectra is plotted in velocity do- 
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TABLE 6 

The line velocities of Ha, H (3, Fen (AA 4924, 5018, 5169) and Hei A5876 as estimated by modelling the observed spectra 

OF SN 2013EJ WITH SYNOW. Fell OR He I LINES VELOCITIES ARE TAKEN TO REPRESENT PHOTOSPHERIC VELOCITY (Uphm). 


UT Date 
(yyyy-mm-dd) 

Phase a 
(day) 

v{ He I ) 

10 3 kms -1 

v( Fen ) 
10 3 kms -1 

v(H«) 

10 3 kms -1 

v(H a ) HV b 
10 3 kms -1 

10 3 kms 1 

v(Hp) HV b 
10 3 kms -1 

2013-08-04.86 

12.1 

8.8 

- 

9.6 

- 

9.7 

- 

2013-08-27.76 

35.0 

- 

6.7 

7.9 

- 

6.6 

- 

2013-09-03.90 

42.1 

- 

5.8 

7.2 

8.5 

5.4 

6.4 

2013-09-29.78 

68.0 

- 

3.6 

5.8 

7.4 

4.0 

5.8 

2013-10-02.89 

71.1 

- 

3.3 

5.4 

7.3 

3.8 

5.8 

2013-10-11.28 

79.5 

- 

3.3 

5.2 

6.3 

3.6 

4.8 

2013-10-27.87 

96.1 

- 

2.7 

4.9 

6.3 

3.5 

4.8 

2013-10-28.79 

97.0 

- 

2.7 

4.8 

6.3 

- 

- 


a With reference to the time of explosion JD 2456497.30 
b High velocity component used to fit the broad Ha and H (3 profile. 



4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 


Rest Wavelength(A) 

Fig. 10.— Comparison of early (12d) and plateau (35d, 68d) 
phase spectra of SN 2013ej with other well-studied type IIP 
SNe 1999em (Leonard et al. 2002a), 1999gi (Leonard et al. 2002b), 
2004et (Sahu et al. 2006; Maguire et al. 2010), 2012aw (Bose et al. 
2013) and 2013ab (Bose et al. 2015). All comparison spectra are 
corrected for extinction and redshift (adopted values are same as 
in Fig. 4). 

main corresponding to rest wavelengths of Hoc At 
12d broad P-Cygni profile (FWHM ~ 9500 kms -1 ) 
is visible which becomes narrower with time as the 
expansion slows down. The blue-shifted absorption 
troughs are direct estimator of expansion velocity of 
the associated line forming layer. The emission peaks 
are found to be blue-shifted (by 3200 kms x at 
12d), which progressively decreases with decrease in 
expansion velocity and almost settling to zero veloc¬ 
ity when the SN starts to enter nebular phase (97d). 
Such blue-shifted emission peaks, especially during early 
phases are generic features observable in SN spec¬ 
tra, e.g., SNe 1987A (Hanuschik & Dachs 1987), 1998A 
(Pastorello et al. 2005), 1999em (Elmhamdi et al. 2003), 
2004et (Sahu et al. 2006), 2012aw (Bose et al. 2013), 



Fig. 11.— The nebular phase spectrum of SN 2013ej at 125d. 
Prominent emission and absorption features are marked and la¬ 
beled. 


2013ab (Bose et al. 2015). These features are tied with 
the density structure of the ejecta, which in turn controls 
the amount of occultation of the receding part of ejecta, 
resulting in biasing of the emission peak (Anderson et al. 
2014a), which are not limited to Ha but applicable to 
all spectral lines. However, such a blue-shift is clearly 
detected for Ha whereas for most other lines, emission 
profiles are weak and peaks are contaminated by adja¬ 
cent P-Cygni profiles. Detailed SN spectral synthesis 
code like cmfgen (Dessart & Hillier 2005b) is capable 
of reproducing such blue-shifted emission peaks. 

As inferred from Fig. 10, the spectral evolution of SN 
2013ej is almost identical to other typical IIP SNe. How¬ 
ever, the comparison of 35 and 68d spectra indicates Fe II 
lines are somewhat under developed as compared to other 
SNe at similar phase. As seen in the 68d comparison, the 
Fell (AA 4924, 5018, 5169) absorption dips are signifi¬ 
cantly weaker in comparison to that seen in other SNe. 

Another prominent and unusual feature is seen in neb¬ 
ular spectra at 109d and 125d, on top of Ha emission, 
and the same is marked as feature A in Fig. 14. This un¬ 
usual dip is resulting into an apparent blue-shift of the 
emission peak, which is in fact larger than that seen in 
the last plateau spectra at 97d. Such evolution is un¬ 
expected and against the general trend of emission peak 
evolution in SNe. The low resolution of these spectra pro¬ 
hibits us from investigating this feature in detail. This 
feature can be split into two emission components, one 
redshifted at 1200 kms -1 and another blueshifted by 
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Fig. 12.— SYNOW modelling of SN 2013ej spectrum at 71d. Model spectrum is shown with thick solid line (blue), while the observed one 
is shown with thin solid line (red). Observed fluxes are corrected for extinction. 




Wavelength(A) Wavelength(A) 

Fig. 13.— SYNOW modelling of SN 2013ej spectra at 8 phases 
during 12d to 97d for H/3, Fell multiplet (left) and Ha (right) 
profiles. Model spectra are shown with thick solid line (blue), while 
the observed ones are shown with thin solid line (red). In the 
model, HI lines are treated as detached to fit the absorption 
troughs. Along with Fell and Hi ; other ions ( Sell , Ball 
, Sill and Nal , Till ) are also incorporated in model to fit 
some weaker features, specially at later phases. In addition to this, 
high-velocity Hi lines are also incorporated (42d onwards) to fit 
the extended Ha and H (3 absorption troughs. The 97d spectrum 
do not have H (3 and Fell wavelength region, hence it is not shown 
here. 

1300 kms -1 (see §A for further explanation) with re¬ 
spect to Ha rest position. Such an asymmetric or double 
peaked Ha nebular emission has been observed in a num¬ 
ber of SNe, e.g. SN 1999em (Leonard et al. 2002a) and 
SN 2004dj (Chugai et al. 2005). Leonard et al. (2002a) 
identified such a dip or notch in Ha emission profile only 
during nebular phase of SN 1999em, which they sug¬ 
gested as possible ejecta-CSM interaction or asymmetry 
in line emitting region. In SN 2004dj, the asymmetry 
in nebular Ha spectra identified by Chugai et al. (2005) 
has been explained by bipolar distribution of 56 Ni with 
a spherical hydrogen envelope (Chugai 2006). 


12. Id 
35. Od 

42. Id 

68. Od 

71. Id 

79.5d 

96. Id 

97. Od 


108.9d 
125.Od 


Fig. 14.— Evolution of Ha line profile at 10 phases during 12d 
to 125d. A zero-velocity line is plotted with a dashed line corre¬ 
sponding to the rest wavelength of Ha A6563. 
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5.4. Ejecta velocity 

Progenitor stars prior to explosion develop stratified 
layers of different elements, which are generally arranged 
in an elemental sequence, hydrogen being abundant in 
the outermost shell, whereas heavier metals like iron pre¬ 
dominate at deeper layers. However at the time of shock 
breakout significant mixing of layers may occur. Spec¬ 
tral lines originating from different layers of the ejecta 
attains different characteristic velocities. Thus study of 
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Fig. 15. — Velocity evolution of Ha, H/3, Hei , Sen and Fen 
lines. The velocities are estimated using blueshift of the absorption 
minima. The expansion velocity of photosphere (n p h m ) estimated 
from SYNOW modeling of He I line at 12d and Fe II lines at later 
phases (see Table 6) are also overplotted for comparison. 

velocity evolution provides important clues to the explo¬ 
sion geometry and the characteristics of various layers. 
Evolution of photospheric layer is of special interest as 
it is directly connected to the kinematics and other re¬ 
lated properties. Photosphere represents the layer of SN 
atmosphere where optical depth attains a value of ~ 2 /% 
(Dessart & Hillier 2005a). Due to complex mixing of lay¬ 
ers and continuous recession of the recombination front, 
no single spectral line can represent the true photospheric 
layer. During the plateau phase, Fell or Sell lines 
are the best estimator of photospheric velocity (u p h). In 
early phases when Fe II lines are not strongly detectable, 
the best proxy for u p h is He I , or H (3 (Takats & Vinko 
2012) in even earlier phases. 

Line velocities can either be estimated by directly lo¬ 
cating the P-Cygni absorption troughs, as done using 
SPLOT task of IRAF, or by by modeling the line profiles 
with velocity as one of the input, as we do in SYNOW. 
In Fig. 15, we plot the line velocities of Ha, H/3, Fell 
(AA 4924, 5018, 5169) and Sen (AA 4670, 6247), using 
the absorption minima method. It is evident that Fell 
and Sc II line velocities are very close to each other and 
they are formed at deeper layers, whereas Ha and H /3 
line velocities are consistently higher at all phases as they 
form at larger radii. The SYNOW estimated photospheric 
velocities are also plotted for comparison, which is very 
close to the Fe II and Sc II velocities estimated from ab¬ 
sorption minima method. Here the SYNOW-derived pho¬ 
tospheric velocities are estimated by modelling He I line 
for 12d spectrum and Fell lines for rest of the spectra. 
Velocities for various lines estimated using SYNOW are 
tabulated in Table 6. 

Fig. 17 shows the comparison of photospheric velocity 
of SN 2013ej with other well-studied type II SNe 1987A, 
1999em, 1999gi, 2004et, 2005cs, 2012aw and 2013ab. For 
the purpose of comparison the absorption trough veloci¬ 
ties have been used, taking the mean of Fe II line triplet, 
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Fig. 16.— For 68d spectrum, the H/3 profile is fitted using synow 
with various velocity components, (a) The fit only with a single 
high velocity component to match the blue wing of the absorption 
dip, (b) with a single low velocity component to match the red 
wing, (c) with single velocity to only fit the trough, (d) with two 
velocity components to fit entire absorption profile. 

or He I lines at early phases where Fe II lines are not de¬ 
tectable. The velocity profile of SN 2013ej is very similar 
to other normal IIP SNe 1999em, 1999gi, 2004et, 2012aw 
and 2013ab, on the other hand velocities of SN 2005cs 
and 1987A are significantly lower. The velocity profile of 
SN 2013ej is almost identical with SNe 2004et, 2012aw 
and 2013ab, whereas it is consistently higher than SNe 
1999gi and 1999em by ~ 800 — 900 km s -1 . For com¬ 
parison of Hi (Ha and H /3) velocities, we have chosen 
all those events which are at least photometrically and 
spectroscopically similar to SN 2013ej. Comparison re¬ 
veals that, H velocities during later phases (60-100 d) 
are consistently higher than all comparable events. SNe 
2012aw and 2013ab, have photospheric velocities identi¬ 
cal to SN 2013ej, but their H velocities are significantly 
lower by large values, e.g., for SN 2013ej the Ha veloc¬ 
ity at 80d is higher by 1500 kms -1 and H /3 is higher by 
2400 kms -1 . Likewise, H velocities for SNe 1999em and 
1999gi are even lower at similar phases. Although SN 
2004et HI velocities are somewhat on higher end, they 
are still significantly less than those of SN 2013ej. It is 
also to be noted that, at 12d SN 2013ej Hi velocities 
are consistent and similar to those of other normal SNe, 
but as it evolves these velocities decline relatively slowly, 
ultimately turning out into a higher velocity profile after 
~ 40d. 

5.5. High velocity components of Hi and CSM 
interaction 

As discussed in §5.2, the broad and extended Ha or H/? 
absorption profiles are not properly reproduced using sin¬ 
gle H i velocity component in SYNOW, and those profiles 
can only be fitted by incorporating a high-velocity (HV) 
components along with the regular one. Fig. 16 shows the 
comparison of SYNOW fits for 68d H /3 profile with various 
single velocity as well as for combined two velocity com- 
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ponents. A single velocity component at 5600 kms -1 
can match the blue wing well and partially the trough, 
whereas, it does not match the red side at all. Simi¬ 
larly, with a single velocity component at 4000 kms -1 
can partially match the red slope of the trough, but does 
not include the trough as well as the extended blue wing. 
By only matching the trough position, the model fits for 
a single velocity of 5300 kms -1 , which does not fit either 
of the blue or red wing. Even-though the ‘detachment’ of 
Hi from photosphere in SYNOW model makes the fit of 
red wing worse by steepening it further, but it is still con¬ 
clusive that none of these single velocity component can 
properly reproduce the absorption profile. It is only by 
including two velocity components together in the model 
could reproduce the entire H/3 profile. Such a scenario 
start to appear from 42d spectrum which only becomes 
stronger as the line evolves until 97d. The Ra troughs 
are also reproduced in a similar fashion. However, it may 
be noted, that such an extended Hi feature may also 
be explained as a possible outcome of a different (com¬ 
plex and extended) density profile which SYNOW can not 
reproduce. 

The comparison of Ra and H/3 velocities with other 
normal IIP SNe (see Fig. 17), estimated by directly lo¬ 
cating the P-Cygni absorption troughs, shows that SN 
2013ej velocities are significantly higher and declines rel¬ 
atively slowly (especially during later phases; 60-100d) as 
compared to those seen in typical IIP SNe, e.g., 1999em, 
1999gi, 2012aw or 2013ab. On the other hand the pho- 
tospheric velocity comparison with other IIP SNe does 
not show any such anomaly. This, we suggest as the ef¬ 
fect of blending with Hi HV components in Ra and 
H/J, which we could separate out while modeling these 
broad features with SYNOW having two velocity compo¬ 
nents. The regular Ra and R/3 velocities estimated from 
SYNOW declines at a normal rate consistent to that seen 
in other SNe (see Fig. 17), whereas the HV components 
remains at higher velocities by 1000 — 2000 kms -1 , de¬ 
clining at relatively slower rate. It is also interesting 
to note that the velocity difference between the regular 
and HV component for Ra and R/3 is similar at same 
epochs. Chugai et al. (2007) identified similar HV ab¬ 
sorption features associated close to Ra and R/3 troughs 
in SNe 1999em and 2004dj, which remained constant 
with time. Presence of such HV features has also been de¬ 
tected in SN 2009bw (Inserra et al. 2012) and SN 2012aw 
(Bose et al. 2013) which is suggestive of interaction of SN 
ejecta with pre-existent CSM. Similar to SN 2013ej, HV 
signatures has been detected all throughout the plateau 
phase evolution of SN 2009bw, while in SN 2012aw such 
features were only detected at late plateau phase (55 to 
104d). Although, we found HV components in SN 2013ej 
by modeling the extended P-Cygni troughs, we are un¬ 
able to visually detect such two individual velocity com¬ 
ponents, this is possibly because of our signal-to-noise- 
ratio limited spectra and weaker strength of HV compo¬ 
nents. Chugai et al. (2007) argued that SN ejecta can in¬ 
teract with the cooler dense shell of CMS material, which 
might have originated from the pre-supernova mass loss 
in the form of stellar winds. Their analysis showed that 
such interaction can led to the detection of HV absorp¬ 
tion features on bluer wings of Balmer lines due to en¬ 
hanced excitation of the outer layers of unshocked ejecta. 
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Fig. 17.— The photospheric velocity (top) evolution (u p h) of SN 
2013ej is compared with other well-studied type II SNe. The u p h 
plotted here are the absorption trough velocities (average of Fell 
lines at late phases and Hel at early phases). Similar comparison 
of P-Cygni absorption velocities, but for Ha and H/3 are shown 
in middle and bottom panels respectively. The regular velocity 
component for Ha and H /3 estimated from SYNOW (without HV 
components; see Table. 6) are also plotted for comparison. 

We, therefore suggest weak or moderate ejecta-CSM in¬ 
teraction in SN 2013ej. X-ray emission from SN 2013ej 
has also been reported by Margutti et al. (2013), which 
they measured a 0.3-10 keV count-rate of 2.7T0.5 cps, 
translating into a flux of ^ 1.1 x 10 -13 erg s -1 cm -2 
(assuming simple power-law spectral model with photon 
index Gamma = 2). Such X-ray emission may also indi¬ 
cate ejecta-CSM interaction suffered by SN 2013ej. 

6. STATUS OF SN 2013ej IN TYPE II DIVERSITY 
6.1. Factors favoring SN 2013ej as type IIL 

Having characterized the event both photometrically 
and spectroscopically, we may now revisit the aspects 
which favor SN 2013ej as type IIL event. The SN 
was originally classified as type IIP (Valenti et al. 2013) 
based on spectroscopic similarity to SN 1999gi. Due to 
same underlying physical mechanisms which govern both 
type IIP and IIL SNe, early spectra may not clearly dis¬ 
tinguish these sub classes of SN type II. The distinguish¬ 
ing factor among IIP and IIL is nominal and mainly de- 
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pend upon light curve characteristics. SN 2013ej shows a 
decline of 1.74 mag (100 d) -1 (see Table 5) or ~ 0.87 mag 
in 50 days, which definitely falls in the criteria of type 
IIL SNe as proposed by Faran et al. (2014). In Fig. 4, 
the spread of template light curves for type IIP and IIL 
(Faran et al. 2014) is shown along with My light curves 
of SNe sample. It is evident that under this scheme of 
classification, SN 2013ej is not a type IIP, rather it is 
marginally within the range of type IIL template light 
curves. This is also justified from the point of basic 
idea behind these classifications, that type IIP must show 
a ‘plateau’ of almost constant brightness for some time 
(^ 90d), which is not the case with SN 2013ej. Due to 
the very fact that SN type II light curves and physical 
properties exhibit a continuum distribution rather than 
a bi-modality (Anderson et al. 2014b), SN 2013ej shows 
intermediate characteristic in the SN type II diversity. 

One distinguishing spectroscopic property Faran et al. 
(2014) found for type IIL SNe is the overall higher pho- 
tospheric ( Fell A5196) velocity and flatter Hi (H/3 and 
Ha) velocity profiles as compared to type IIP counter¬ 
part. Although Fell velocities are on the higher end as 
compared to typical IIP SNe velocities, we do not find it 
as a remarkable deviation to distinguish SN 2013ej from 
IIP sample. However, we do see a anomaly in Ha, H /3 
absorption minima velocity profiles, as they start off with 
velocities consistent with those of type IIP but declines 
relatively slowly (see §5.4 for more description of this fea¬ 
ture) ultimately surpassing faster declining IIP velocity 
profiles after 50 days. This characteristic feature of HI 
velocities for SN 2013ej is typical for most IIL SNe as 
found by Faran et al. (2014). 

6.2. CSM interaction and type IIL 

Faran et al. (2014) proposed a possible explanation for 
the hatter velocity profiles in IIL SNe, which is due the 
lack of hydrogen in deeper and slow expanding layers 
of ejecta, resulting into higher Hi absorption veloci¬ 
ties arising mostly from outer layer. However, for SN 
2013ej we suggest the battening of Ha and H f3 veloc¬ 
ity profiles are due to the contamination of HV compo¬ 
nent of Hi (see §5.5). Indication of CSM interaction 
in SN 2013ej may also be inferred from X-ray detection 
by Margutti et al. (2013). Valenti et al. (2015) found SN 
2013by, a type IIL SN, to be moderately interacting with 
CSM. This led them to ask the prevalence of CSM inter¬ 
action among IIL SNe in general. Type IIL SNe orig¬ 
inate from progenitors similar to IIPs, but have lost a 
significant fraction of hydrogen before explosion during 
pre SN evolution. Hence it may not be usual to detect 
HV Hi signatures in Ha, H f3 absorption profiles as a 
consequence of ejecta-CSM interaction. A moderate or 
weak interaction may produce a HV component blend¬ 
ing with Ha, H/3 profiles, which may result into shift in 
absorption minima, rather than a prominent secondary 
HV dip. Such a scenario may perfectly explain the rel¬ 
atively higher and hatter Hi velocity prohles of most 
type IIL SNe as compared to IIP counterparts, found by 
Faran et al. (2014) based on direct velocity estimates of 
absorption minima. 

Another example of CSM interaction in type IIL is SN 
2008fq, which does show strong interaction signature like 
a type Iln (Taddia et al. 2013), but also shows a steep 
decline like IIL during hrst 60 days (Faran et al. 2014). 


Supernova PTFlliqb (Smith et al. 2015) is also a type 
Iln SN, having prominent CSM interaction signatures, 
but with IIL like steeper light curve. Initial spectra of 
this SN showed Iln characteristics, however late plateau 
spectra revealed features similar to type IIL. PTFlliqb 
originated from a progenitor identical to type IIP/L, in¬ 
stead of a LBV as expected for a typical Iln. However, 
because of rare detection of type IIL events and its fast 
decline in magnitudes we do not have sufficient informa¬ 
tion to investigate CSM interaction in ah such objects. 
Thus, the question still remains open if ah or most IIL 
SNe interact with CSM and whether the hatter HI ab¬ 
sorption minima velocity prohles is a consequence of in¬ 
teraction. 


7. LIGHT CURVE MODELLING 

To determine the explosion parameters of SN 
2013ej, the observed light curve is modeled follow¬ 
ing the semi-analytical approach originally developed 
by Arnett (1980) and further rehned in Arnett & Fu 
(1989). More appropriate and accurate approach 
would have been detailed hydrodynamical modeling (e.g. 
Falk & Arnett 1977; Utrobin 2007; Bersten et al. 2011; 
Pumo & Zampieri 2011) to determine explosion proper¬ 
ties, however application of simple semi-analytical mod¬ 
els (Arnett 1980, 1982; Arnett & Fu 1989; Popov 1993; 
Zampieri et al. 2003; Chatzopoulos et al. 2012) can be 
useful to get preliminary yet reliable estimates of the 
parameters without running resource intensive and time 
consuming hydrodynamical codes. Nagy et al. (2014) 
also followed the original semi-analytical formulation pre¬ 
sented by Arnett & Fu (1989) and modeled a few well 
studied II SNe. The results are compared with hydrody¬ 
namical models from the literature and are found to be in 
good agreement. The model light-curve is computed by 
solving the energy balance of the spherically symmetric 
supernova envelope, which is assumed to be in homolo¬ 
gous expansion having spatially uniform density profile. 

The temperature evolution is given as (Arnett 1980), 

T(x,t) 4 = T$ip(x)(p{t) > 

where x is defined as dimensionless co-moving radius rel¬ 
ative to the mass of the envelope and, 'ip(x) is the radial 
component of temperature profile which falls off with ra¬ 
dius as sin(7rx)/7rx. Here we incorporate the effect of 
recombination, as shock heated and ionized envelope ex¬ 
pands and cools down to recombine at temperature T rec . 
We define Xi as the co-moving radius of the recombi¬ 
nation front and the opacity (/c) changes very sharply 
at this layer such that n ~ 0 for the ejecta above X{. 
Following the treatment of Arnett & Fu (1989) the tem¬ 
poral component of temperature, c/)(t) can be expressed 
as (Nagy et al. 2014), 


d<j>(t) 

dz 


m 

Roxf 


PlC(t) ~P2<l>(t)Xi - 2 x?<p(t) 


Ro dxi 


R(t) dz 


here £(£) is the total radioactive energy input from de¬ 
cay chain of unit mass of 56 Ni, which is normalized to 
the energy production rate of 56 Ni. The rest of the pa¬ 
rameters in the equation have usual meaning and can 
be found in aforementioned papers. From this ordinary 
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Fig. 18. — Model fit (solid line) on the observed bolometric light 
curve (open circles) of SN 2013ej. The green solid line follows 
only the radioactive decay law, where the recombination front has 
completely disappeared. 

differential equation we can find out the solution of (/>{t) 
using Runge-Kutta method. The treatment adopted to 
determine Xi is somewhat similar to Nagy et al. (2014), 
where we numerically determine the radius Xi (to an ac¬ 
curacy of 10 -12 ) for which the temperature of the layer 
reaches T rec . Once we find out the solution of </>(£) and 
x*, the total bolometric luminosity is calculated as the 
sum of radioactive heating and rate of energy released 
due to recombination, 

L(t) = x M^ Eth ^ (i _ e~ Ag/t2 ^\+4nri Qp(xi,t)R(t)^-, 

here, d(xi)/dt is the inward velocity of co-moving recom¬ 
bination front and the term [1 — exp(—A g /t 2 )\, takes into 
account of gamma-ray leakage from the ejecta. The fac¬ 
tor A g is the effectiveness of gamma ray trapping (see 
e.g., Clocchiatti & Wheeler 1997; Chatzopoulos et al. 
2012), where large A g means full trapping of gamma 
rays, this factor is particularly important to model the 
SN 2013ej tail light curve. In this relation we also modi¬ 
fied the second term to correctly account for the amount 
of envelope mass being recombined. 

To model SN light curves it is essential to obtain the 
true bolometric luminosity from observations. Since our 
data is limited only to optical and UV bands, we adopt 
the prescription for color dependent bolometric correc¬ 
tions by Bersten & Hamuy (2009) to obtain bolometric 
light curve for SN 2013ej. Figure 18 shows the model 
fit with the observed bolometric light curve of the SN. 
We estimate an ejecta mass of 12 M 0 , progenitor radius 
of 450 R 0 and explosion energy (kinetic ± thermal) of 
2.3 foe (10 51 erg). The uncertainty in mass and radius is 
about 25%. We find that the plateau duration is strongly 
correlated with explosion energies (especially kinetic), 
and also with k and T rec . Thus depending upon these 
parameters our model is consistent with a wide range of 
explosion energies, with 2.3 foe towards the lower end 
and energies up to 4.5 foe at higher end. Assuming the 
mass of the compact remnant to be 1.5-2.0 M 0 , the total 
progenitor mass adds up to be 14M 0 . 

The mass of radioactive 56 Ni estimated from the model 


is O.O18M 0 , which primarily governs the tail light curve 
of the SN. As discussed in §4.2, the slope of the tail light 
curve observed for SN 2013ej is significantly higher than 
other typical IIP SNe and also to that expected from ra¬ 
dioactive decay of 56 Co to 56 Fe. The light curve powered 
by full gamma-ray trapping from radioactive decay chain 
of 56 Ni —>> 56 Co —» 56 Fe results in a slower decline and 
does not explain the steeper tail observed in SN 2013ej. 
In the model we decreased the gamma-ray trapping effec¬ 
tiveness parameter A g to 3 x 10 4 day 2 , which matches the 
steeper radioactive tail. The gamma-ray optical depth 
can be related to this parameter as r g ~ A g /t 2 . This 
implies that the gamma-ray leakage in SN 2013ej is sig¬ 
nificantly higher than other typical type IIP SNe. 

Valenti et al. (2014) using early temperatures (< 5 
days) of SN 2013ej provided a preliminary estimate of the 
progenitor radius as 400 — 600 R 0 , which is in good agree¬ 
ment with our result. Our progenitor mass estimate is 
also consistent with that reported by Fraser et al. (2014) 
from direct observational identification of the progenitor 
using HST archival images, which is 8 — 15.5 M 0 . 

8. SUMMARY 

We present photometric and spectroscopic observa¬ 
tions of SN 2013ej. Despite low cadence optical photo¬ 
metric follow up during photospheric phase, we are able 
to cover most of the important phases and features of 
light curve. 

Our high resolution spectrum at 80d shows the pres¬ 
ence of Nal D (AA 5890, 5896) doublet for Milky Way, 
while no impression for host galaxy NGC 0628. This in¬ 
dicates that SN 2013ej suffers minimal or no reddening 
due to its host galaxy. 

The optical light curves are similar to type IIL SNe, 
with a relatively short plateau duration of 85d and 
steeper decline rates of 6.60, 3.57, 1.74, 1.07 and 0.74 mag 
100 day -1 in UBVRI bands respectively. The compari¬ 
son of absolute V band light curves shows that SN 2013ej 
suffers the higher decline rate than all type IIP SNe, but 
similar to type IIL SNe 1980k, 2000dc and 2013by. The 
drop in luminosity during the plateau-nebular transition 
is also higher than most type II SNe in our sample, which 
is 2.4 mag in V band. 

The UVOT UV optical light curves shows steep decline 
during first 30 days at a rate of 0.182, 0.213, 0.262 mag 
d -1 in uvwl, uvw2 and uvm2 bands respectively. The 
absolute UV light curves are identical to SN 2012aw and 
also shows a similar UV-plateau trend as observed in SN 
2012aw. 

Owing to the large drop in luminosity during plateau- 
nebular transition, the light curve settles to a signifi¬ 
cantly low luminous tail phase as compared to other 
normal IIP SNe. The mass of radioactive 56 Ni esti¬ 
mated from the tail bolometric luminosity is 0.020±0.002 
M 0 , which is in between normal IIP SNe (e.g., 1999em, 
2004et, 2012aw) and subluminous events, like SN 2005cs. 

The spectroscopic features and their evolution is simi¬ 
lar to normal type II events. Detailed SYNOW modelling 
has been performed to identify spectral features and to 
estimate velocities for Hcq H/?, Fell (AA 4924, 5018, 
5169) and Sell (AA 4670, 6247) lines. The photospheric 
velocity profile of SN 2013ej, which is represented by Fe II 
lines and Hel line at 12d, is almost identical to SNe 
2004et, 2012aw and 2013ab. The Hcq H/3 velocities es- 
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timated by directly locating the absorption troughs are 
significantly higher and slow declining as compared to 
other normal IIP events. However, such Hi velocity 
profiles are typical for type IIL SNe. 

The P-Cygni absorption troughs of Ha and H /3 are 
found to be broad and extended which a single HI com¬ 
ponent in SYNOW model could not fit properly. However, 
these extended features are fitted well with SYNOW by in¬ 
corporating a high velocity Hi component. These HV 
components can be traced throughout the photospheric 
phase which may indicate possible ejecta-CSM interac¬ 
tion. Our inference is also supported by the detection of 
X-ray emission from the SN 2013ej (Margutti et al. 2013) 
indicating possible CSM interaction, and the unusually 
high polarization reported by Leonard et al. (2013) may 
also further indicate asymmetry in environment or ejecta 
of the SN. Such CSM interaction and their signature in 
Ha, H ft profiles has also been reported for SNe 2009bw 
(Inserra et al. 2012) and 2012aw (Bose et al. 2013). 

Nebular phase spectra during 109 to 125d phases are 
dominated by characteristic emission lines, however the 
Ha line shows an unusual notch, which may be explained 
by superposition of HV emission on regular Ha profile. 
Although, the origin of the feature is not fully explained, 
it may indicate bipolar distribution of 56 Ni in the core. 

We modeled the bolometric light curve of SN 2013ej 
and estimated a progenitor mass of ~ 14M 0 , radius of 
~ 45OR 0 and explosion energy of ^ 2.3 foe. These pro¬ 
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Fig. 19.— Ha profile of 125d spectrum fitted by two component gausian profile seperated by ~2500 kms 1 . 


APPENDIX 

NEBULAR Ha PROFILE 

The unusual notch in nebular Ha profile can be described as a superimposition of two profiles. Fig. 19 shows the 
observed Hct profile at 125d which is fitted by two component Gaussian profiles. These two profiles are separated by 
55 A 2500 kms 1 ), one being blue shifted by -1300 kms 1 while the other is red shifted at 1200 kms 1 with 
respect to rest Ha position. The FWHM for the blue component is 54A and for red component is 146A. The redshifted 
component is dominant in strength over the blue one, having their ratio of equivalent widths to be 4.5. It may be 
noted that for the sake of simplicity and only for the purpose of illustration we used Gaussian profiles, which does not 
account for the P-Cygni absorption troughs as we see on bluer wings of line profiles in observed SN spectrum. 











